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FORECASTING  THE  RETENTION  OF  NAVY  PILOTS 
by 

S . Zacks 


1.  Introduction 


The  present  paper  addresses  the  problem  of  estimating  the  number  of 
Navy  pilots  (or  other  career  force  officers  of  a certain  profession)  that 
are  expected  to  remain  in  service.  More  specifically,  we  consider  in  the 
present  study  the  retention  of  Navy  pilots  who  are  serving  their  first  term 
in  an  operational  squadron,  after  completing  training  schools.  Generally, 
during  the  first  term  of  operational  service  there  is  very  little  attrition. 
Attrition  at  this  time  is  due  mainly  to  accidents  or  sickness.  After  the 
first  term  of  service  the  withdrawal  of  pilots  from  operational  service 
becomes  quite  significant.  It  is  highly  Important  to  develop  a method  for 
forecasting  how  many  pilots  from  each  cohort  will  remain  in  service.  As 
will  be  shown  later,  the  retention  of  pilots  in  service  depends  to  a large 
extent  on  their  original  service.  These  sources  are  coded  by  certain 
original  source  codes  (OSC).  In  our  study  we  analyzed  data  from  six  such 
original  sources  (OSC  1,  OSC  3,  OSC  4,  OSC  6,  OSC  38  and  OSC  80). 

In  addition,  the  data  was  classified  according  to  the  date  of  the 
minimum  service  requirement  initials  (MSRI)  of  the  pilots  and  how  many 
years  they  were  before  or  after  their  MSRI.  In  Section  2 we  provide  six 
tables  of  the  actual  data.  In  Section  3 we  present  the  statistical  model 
and  the  methodology  of  forecasting  the  retention  of  the  officers.  This 


methodology  is  based  on  determining  tolerance  limits  for  estimating  predic- 
tion intervals  of  the  anticipated  number  of  personnel  retained  in  the  service. 
In  Section  4 we  present  the  numerical  results.  A Fortran  computer  program 
is  given  in  the  appendix.  We  would  like  to  comment  at  this  point  that  the 
forecasting  procedure  developed  here  is  not  restricted  only  to  Navy  pilots. 

It  is  expected  to  perform  well  if  applied  to  other  professional  groups  of 
career  force  officers.  The  theoretical  development  in  the  present  paper  is 
similar  to  that  applied  in  the  forecasting  of  the  total  size  of  the  Marine 
Corps  (see  Haber  and  Zacks  [1],  Zacks  [2],  and  Zacks  [3]). 

2.  The  Data 

The  data  consist  of  six  tables  (Table  1 - Table  6)  in  which  the  number 
of  pilots  in  service  is  tabulated  according  to  the  quarter  of  the  year  they 
attain  MSRI  and  the  time  period  (in  years)  before  or  after  the  MSRI.  Each 
table  corresponds  to  a different  OSC,  and  spans  the  period  from  the  begin- 
ning of  1971  to  the  middle  of  1975.  For  the  computer  application  we  opened 
six  data  files  (PIL  1,  PIL  3,  PIL  4,  PIL  6,  PI  38  and  PI  80)  consisting  of 
eighteen  lines  (card).  Each  line  (card)  contains  eight  three-digit  numbers. 
Lines  110  - 160  and  230  - 310  of  the  programs  in  the  appendix  provide  the 
instructions  for  calling  these  files  and  reading  the  data  in  the  specified 
manner. 

3.  The  Forecasting  Procedure 

3.1  The  Statistical  Model 

Let  X(l,J,t)  denote  the  number  of  pilots  from  the  1th  source 
(1  ■ 1,...,6)  , with  specified  MSRI  at  the  jth  quarter  (J  ■ 1,...,18) 
who  are  still  in  service  at  time  t , t - 1,...,8  . The  ...  MSRI  quarters 
are  specified  in  the  following  table. 
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The  time  index  t corresponds  to  the  number  of  years  before  (-)  or  after 
(+)  the  MSRI  (o)  in  the  following  manner: 

tl2345678 
Years  -2 . E3  0 O T.  TTs  IT 

For  a specific  combination  (i,j)  , given  the  value  of  X(i,j,t)  we  wish  to 

predict  the  values  that  this  variable  will  assume  at  time  t+v  , v = If 2, 

3,4,...  . The  statistical  model  assumes  that,  given  the  value  of  X(i,j,t)  , 

the  conditional  distribution  of  X(i,j,  t+v)  is  binomial  over  the  range 

{0,...,  X(i,j,t)}  . More  specifically,  given  the  observed  value  of  X(i,j,t)  , 

n^j^  say,  the  future  value  of  X(i,j,  t+v)  , is  a random  variable  having  a 

(conditional)  binomial  distribution,  B(n. . , 9.^  ) , where  6,^  is  the 

ijt  itv  itv 

retention  probability  for  the  period  (t,  t+v)  . Generally,  we  say  that  a 
random  variable  X has  a binomial  distribution  B(N,d)  if  its  probability 
function  is 

P[X  - JJ  « (j)  (1-0)”"^  , j - 0,1 N . (3.1) 

We  assume  here  that  these  retention  probabilities  do  not  depend  on  the  MSRI 
quarter  of  Che  cohort  under  consideration.  In  other  words,  all  cohorts  of 
pilots  coming  from  the  same  OSC  at  different  time  periods  have,  at  similar 
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TABLE  1 


NUMBER  OF  NAVY  PILOTS  FROM 
OSC  1 IN  SERVICE  (FILE  PIL  1) 


MSRI 

Year 

Date 

Quarter 

1 

f 

Nuniber 

of  Years  Before  or 

After 

MSRI 

- - 

-1. 

-.5 

0 

0.5 

1. 

1.5 

2. 

1971 

1 

i 

( 

1 57 

57 

57 

48 

42 

42 

41 

41 

2 

25 

25 

24 

20 

19 

19 

18 

16 

3 

' 22 

21 

21 

15 

14 

14 

13 

12 

4 

97 

97 

96 

85 

78 

72 

68 

58 

1972 

1 

47 

47 

46 

41 

41 

37 

33 

30  1 

22 

22 

22 

21 

18 

17 

15 

14  1 

57 

56 

56 

53 

48 

45 

37 

36  j 

125 

124 

122 

111 

102 

95 

80 

78  1 

1973 

1 

57 

55 

55 

51 

43 

35 

35 

r 

i 

1 

2 

24 

23 

23 

22 

21 

20 

17 

1 

t 

3 

8 

8 

8 

7 

7 

7 

i 

4 

177 

176 

171 

155 

130 

124 

1 

1974 

1 

33 

31 

30 

27 

27 

1 

i 

2 

12 

12 

12 

10 

8 

3 

9 

9 

8 

8 

1 

4 

74 

74 

72 

62 

f 

» 

1975 

1 

66 

66 

65 

2 

56 

54 

54 

‘ 
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TABLE  2 


NUMBER  OF  NAVY  PILOTS  FROM 
OSC  3 IN  SERVICE  (FILE  PIL  3) 


MSRI 

Year 

Date 

Quarter 

Number 

of  Years  Before  or 

After 

MSRI 

-2. 

-1. 

-.5 

0 

0.5 

1. 

1.5 

2. 

1971 

1 

u 

118 

071 

047 

046 

044 

039 

038 

2 

!140 

085 

077 

060 

058 

052 

050 

042 

3 

|l71 

140 

137 

084 

073 

066 

058 

053 

4 

*117 

113 

107 

063 

060 

051 

044 

042 

1972 

1 

192 

189 

182 

120 

110 

089 

084 

077 

2 

258 

244 

224 

144 

128 

108 

097 

087 

3 

243 

234 

186 

122 

105 

091 

076 

070 

4 

179 

169 

150 

091 

076 

070 

061 

059 

1973 

1 

181 

146 

130 

083 

066 

057 

049 

2 

296 

243 

207 

137 

121 

108 

105 

3 

162 

125 

113 

077 

065 

064 

4 

227 

172 

145 

111 

092 

088 

1974 

1 

235 

173 

160 

109 

094 

2 

341 

242 

240 

184 

177 

3 

171 

132 

127 

090 

4 

171 

152 

134 

116 

1975 

1 

179 

172 

167 

2 

174 

173 

170 

- 5 - 
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TABLE  3 


NUMBER  OF  NAVY  PILOTS  FROM 
OSC  4 IN  SERVICE  (FILE  PIL  4) 


MSRI 

Year 

— 

Date 

Quarter 

Numb e I 

oC  Years  Before  or 

After 

MSRI 

-2. 

-1. 

-.5 

0 

0,5 

1. 

1.5 

2. 

1971 

1 

037 

036 

035 

028 

022 

021 

020 

019 

2 

035 

035 

032 

025 

024 

023 

023 

020 

3 

066 

066 

064 

051 

046 

046 

042 

040 

1 

4 

095 

089 

088 

071 

065 

061 

057 

050 

1972 

1 

042 

040 

040 

034 

029 

024 

023 

020  1 

2 

038 

037 

037 

034 

032 

031 

028 

026  1 

3 

056 

056 

056 

052 

049 

046 

034 

034  1 

4 

078 

074 

074 

065 

058 

051 

041 

036  1 

1973 

1 

041 

040 

040 

032 

031 

023 

018 

1 

2 '047 

047 

045 

036 

031 

026 

025 

3 

093 

092 

088 

07^ 

057 

052 

4 

037 

037 

035 

032 

024 

021 

1974 

1 

012 

012 

012 

Oil 

010 

. 

2 

024 

024 

024 

018 

015 

1 

3 

064 

059 

057 

043 

i 

4 

060 

055 

052 

046 

I 

1975 

1 

032 

031 

028 

1 

2 

034 

032 

032 

V . 


J 
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TABLE  4 


NUMBER  OF  NAVY  PILOTS  FROM 
OSC  6 IN  SERVICE  (FILE  PIL  6) 


MSRI 

Year 

Date 

Quarter 

Number  of  Years  Before  or 

After  MSRI 

1 

-2. 

-1. 

-.5 

0 

0.5 

1. 

1.5 

2.  j 

1971 

1 

Oil 

Oil 

007 

006 

005 

003 

003 

003  ^ 

2 

021 

019 

018 

013 

012 

Oil 

010 

009  ' 

3 

022 

022 

021 

013 

013 

013 

013 

013 

4 

016 

016 

015 

Oil 

010 

009 

009 

009  i 

1972 

1 

016 

016 

015 

008 

007 

00/ 

005 

005 

2 

020 

017 

017 

012 

Oil 

OC/} 

008 

007  I 

3 

017 

015 

013 

010 

010 

0oc5 

008 

008  1 

4 

007 

007 

006 

004 

003 

003 

003 

003  j 

1973 

1 

018 

015 

013 

008 

007 

007 

007 

t 

2 

033 

029 

027 

018 

014 

013 

013 

1 

3 

022 

015 

014 

010 

010 

009 

4 

008 

005 

004 

003 

002 

001 

J 

1974 

1 

006 

004 

004 

002 

1 

2 

006 

003 

003 

002 

002 

j 

3 

005 

004 

003 

002 

4 

006 

006 

005 

005 

1975 

1 

002 

002 

002 

2 

002 

002 

002 

i 
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TABLE  5 

NUMBER  OF  NAVY  PILOTS  FROM 
OSC  38  IN  SERVICE  (FILE  PL  38) 


MSRI 

Year 

Date 

Quarter 

i 

Number 

of  Years  Before  or 

After 

MSRI 

i-2. 

-1. 

-.5 

0 

0.5 

1. 

1.5 

2. 

1971 

1 

014 

014 

010 

007 

008 

008 

008 

008 

2 

023 

016 

016 

Oil 

Oil 

010 

010 

009 

3 

014 

014 

014 

Oil 

010 

010 

008 

008 

4 

006 

006 

005 

003 

003 

002 

002 

002 

1972 

1 

005 

005 

005 

002 

001 

001 

001 

2 

017 

017 

017 

013 

013 

Oil 

Oil 

010 

3 

013 

012 

Oil 

008 

007 

007 

006 

006 

4 

004 

003 

003 

002 

002 

002 

002 

002 

1973 

1 

020 

018 

016 

012 

Oil 

010 

010 

1 

2 

017 

016 

013 

007 

007 

007 

007 

3 

013 

Oil 

010 

009 

008 

008 

4 

012 

010 

009 

006 

006 

006 

1974 

1 

013 

010 

010 

008 

007 

2 

021 

018 

018 

013 

013 

3 

008 

007 

006 

006 

4 

007 

006 

004 

003 

1975 

1 

002 

001 

001 

2 

001 

001 

001 

i 

- 8 - 
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TABLE  6 


NUMBER  OF  NAVY  PILOTS  FROM 
OSC  80  IN  SERVICE  (FILE  PL  80) 


MSRI 

Year 

Date 

Quarter 

I 

1 

Number 

of  Years  Before  or 

After 

MSRI 

' -2. 

-1. 

-.5 

0 

0.5 

1. 

1.5 

2. 

1971 

1 

Rj 

■1 

001 

001 

2 

002 

002 

mm 

001 

001 

3 

,005 

PIOCl 

4 

Oil 

008 

006 

006 

006 

005 

004  j 

1972 

1 

006 

005 

004 

002 

002 

001  j 

2 

005 

003 

003 

003 

003 

003  . 

3 

019 

013 

017 

015 

013 

013 

010 

010  1 

4 

040 

036 

031 

026 

023 

019 

016 

013  i 

1973 

1 

022 

020 

018 

013 

010 

010 

008 

1 

i 

2 

021 

019 

019 

016 

014 

013 

013 

3 

041 

036 

024 

018 

018 

1 

j 

4 

064 

045 

036 

030 

024 

022 

i 

i 

1974 

1 

034 

wM 

027 

019 

019 

2 

WBBM 

026 

021 

021 

t 

3 

018 

017 

014 

4 

016 

013 

Oil 

i 

1975 

1 

031 

028 

027 

2 

055 

054 

055 

T-358 


time  points,  t , of  tlieir  service,  the  same  retention  probabilities. 
Moreover,  we  assume  that,  for  eacli  i = , and  each  t = ; 

{X(i,J,t)  ; j “ 1,...,18}  are  independent  random  variables. 


3.2  Confidence  Limits  For  Retention  Probabilities 


It  is  well  known  that  the  sum  of  independent  binomial  random  variables, 
with  the  same  probability  of  success  p , has  a binomial  distribution  with 
the  same  probability  of  success  p and  number  of  trials  which  is  the  sum  of 
the  corresponding  ones.  Accordingly,  the  conditional  distribution  of  the 
18 

sun  T(l,  t+v)  = I X(i,j,  ttv)  given  {X(i,l,t)  “ n^j^^ , . . . ,X(i,18,t)  - 

18 

is  the  binomial  distribution  B(N.  ,6,  ) , where  N,  = 2 ,, 

it  itv  it  n^j^  . We 

estimate  by  determining  (1-a)  confidence  limits  for  this  parameter 

based  on  the  observed  values  of  T(l,  t+v)  and  N . ^ . Let  9,^  = 

’ it  itv 

(T(i,  t+v)  + -5)/N^^  + 1)  . For  large  values  of  , the  transformed 


variable  Y.  =2  arcsin 
itv 


6^^^  is  approximately  normally  distributed  with 


mean  n 


itv 


= 2 arcsin  9 


itv 


and  variance  D 


it 


1/N 


it 


(See  Johnson  and 


Kotz  t^])»  Accordingly 


(1) 

/ / 

1 — \ 

. « max 

itv 

^itv  " ^l-a/2/ 

(2) 
itv  ’ 

(’'•  ’'^Itv  ^1-0/2/ 

(3.2) 


arc  lower  and  upper  (l-oi)  confidence  limits  for  the  transformed  parameter 
Pltv  ; where  bhe  (l-a/2)  fractlle  of  the  standard  normal  dis- 

tribution. Furthermore,  since  the  2 arcsin  17  transformation  is  strictly 
increasing  over  the  range  of  0 ^ 9 ^ 1 , (1-a)  confidence  limits  for  the 
retention  probability  obtained  from  (3.2)  by 
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(k) 


0 


Itv 


k = 1,2  . 


(3.3) 


(1)  (2) 

and  are  the  lower  and  upper  confidence  limits,  respectively. 


3. 3 Prediction  Intervals  For  Forecasting 

On  the  basis  of  these  confidence  limits  we  compute  a prediction  in- 
terval for  X(l,j,  t+v)  given  that  X(i,j,t)  « *^ijt  ‘ ^ mentioned  earlier, 

the  conditional  distribution  of  X(i,j,  t+v)  given  that  X(i,j,t)  = n^j^. 

is  the  binomial  » ®itv^  ' ®itv  known,  we  anticipate 

that  with  probability  y the  random  variable  X(l,j,  t+v)  will  assume 
values  in  the  interval 


The  limits  of  this  interval  are  defined  by 


BU(p;N,0)  = least  non-negative  integer, 

(3.5) 

j , such  that  Pq[X  1 j]  ^ p ; 
where  X has  the  B(N,6)  distribution;  and 


BL(p;N,0)  “ maximal  non-negative  integer. 

J , such  that  P„ [X  < 1 ] < p . 

" y ” 


(3.6) 


The  interval  (3.4)  is  called  a y - content  prediction  interval.  In  almost 
all  practical  cases  the  retention  probabilities  unknown.  We 

therefore  estimate  the  prediction  Interval  (3.4)  by  substituting  the  lower 

(1) 

confidence  limit  0.^  for  0.  In  (3.6)  and  the  upper  confidence  limit 

(2) 

0,  for  0.  In  (3.3).  The  resulting  Interval  contains,  with  confidence 
Itv  Itv 
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probability  1-a  , the  true  prediction  interval  (3.4).  It  is  called  a 
(l-a,Y)  tolerance  interval. 


The  determination  of  the  tolerance  limits  BL 

(2)\ 


in. 


(1) 

”ijt  ’ ®itv 


) 


and  BU 


itL  . 


"ijt  ’ °itVj 


follows  an  algorithm  based  on  the  evaluation 


of  the  binomial  cumulative  distribution  functions 

(2)\ 


and  B 


fijt  * 


(c.d.f.)  B^n^j^  . 


itv  i 


These  c.d.f. 's  are  computed  by  a subroutine  FUNCTION 


CB(I,N,P)  given  in  program  NAVYS  of  the  appendix.  If  N is  greater  than 
50  the  binomial  c.d.f.  CB(I,N,P)  is  approximated  by  a normal  c.d.f.,  according 
to  the  well-known  approximation 


CB(j,N,p)  4, 


,1+.5  - Np 

Jnp(I-p) 


(3.7) 


where  the  function  ^(O  is  the  standard  normal  c.d.f..  This  standard  normal 
c.d.f.  is  determined  by  the  subroutine  FUNCTION  CNDX(Y)  (in  lines  1820  - 2010 
of  program  NAVYS).  The  algorithm  of  computation  is  constructed  according  to 
Formulae  (3.5)  and  (3.6)  (see  subroutine  function  KBL(P,N,T)  and  KBU(P,N,T). 

The  program  is,  however,  relatively  slow  since  the  subroutine  functions  require 
many  repetitive  computations.  We  have,  therefore,  applied  the  normal  approxima- 
tion to  the  binomial  distribution  even  if  N is  not  large,  and  applied  the 
approximation 


Program  NAVYR  is  based  on  this  approximation.  It  turned  out.  Indeed,  that 
program  NAVYR  determines  the  tolerance  limits  about  five  times  faster  than 
program  NAVYS.  The  tolerance  limits,  however,  do  not  differ  in  most  of  the 
cases  in  more  than  one  unit.  In  Table  7 we  provide  the  tolerance  limits 
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Table  7 

Tolerance  Limits  (y“.95,  l-a=.95) 
by  Programs  NAVYS  and  NAVYR:  OSC  1 


NAVYS 

NAVYR 

j/t 

■1 

5 

6 

7 

8 

8 

7 

6 

5 

4 

! 9 

26 

26 

t • 

11 

11 

10 

14 

13 

3 

4 

11 

7 

6 

6 

6 

0 

0 

0 

1 

1 

12 

89 

83 

83 

89 

i 

57 

50 

50 

58 

13 

24 

22 

21 

20 

22 

24 

i 

13 

9 

7 

8 

9 

13 

14 

9 

8 

8 

8 

8 

9 

1 

i 

2 

1 

1 

1 

2 

3 

15 

8 

8 

8 

7 

7 

7 

8 

8 

( 

1 

2 

1 

0 

1 

1 

2 

2 

16 

49 

53 

47 

44 

43 

46 

53 

48 

i 

1 

29 

34 

26 

22 

23 

26 

35 

30 

i 

17 

58 

51 

55 

49 

46 

45 

48 

55 

51 

58  ! 

40 

31 

36 

27 

23 

24 

28 

37 

32 

41  j 

18 

49 

43 

47 

41 

39 

38 

41 

46 

43 

48  i 

» 

33 

25 

29 

22 

19 

19 

22 

30 

26 

33 

1 
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determined  by  programs  NAVYS  and  NAVYR  for  the  data  of  Table  1.  As  seen 
In  Table  7,  the  approximation  provided  by  program  NAVYR  is  quite  good  and 
justifies  the  use  of  this  program,  which  saves  a considerable  amount  of 
computing  time. 

4.  Numerical  Forecasting  of  Pilot  Retention 

In  the  present  section  we  provide  the  numerical  results  obtained  by 
applying  program  NAVYR  to  the  data  of  Tables  1-6.  The  results  are  given 
in  Tables  8-11.  Each  table  extends  over  ten  values  of  j (j  " 9,..., 18) 
which  corresponds  to  cohorts  entering  from  the  first  quarter  of  1973.  There 
are  eight  columns  corresponding  to  values  of  t from  2 (one  year  before 
MSRI)  to  8 (2  years  after  MSRI).  For  each  value  of  j we  present  the  upper 
tolerance  limit  (first  row),  lower  tolerance  limit  (second  row),  and  actual 
value  (third  row).  In  all  cases  where  the  actual  values  have  not  yet  been 
observed  we  print  the  value  o . These  cases  appear  on  the  lower  right  side 
of  each  table.  The  other  cases  correspond  to  past  periods.  The  tolerance 
limits  were  computed  for  past  periods  for  the  sake  of  testing  whether  the 
actual  values  fall  in  the  tolerance  intervals.  As  seen  in  the  tables,  the 
tolerance  Intervals  indeed  cover  most  of  the  actual  values.  Finally,  the 
program  adds  the  upper  and  lower  tolerance  limits  of  all  cohorts  for  one, 
two,  three,  and  four  periods  ahead  to  provide  the  prediction  Intervals  for 
the  number  of  pilots  from  those  presently  in  service,  from  all  the  six  OSC's, 
which  will  be  in  service  6,  12,  18,  and  24  months  from  the  forecasting  date. 
For  the  present  data  under  consideration  these  limits  are 

Table  12 

Tolerance  Limits  For  All  Pilots  From  The  6 OSC's 


6m 

12m 

18m 

24m 

Upper 

1519 

1257 

988 

686 

Lower 

822 

662 

537 



362 

14  - 


T-358 


Table  9 

Tolerance  Limits  For  OSC  3 Cohorts 
(r=.95,  l-a=.95) 


j/t 


2 


3 


4 


11 

51 

41 

29 

37 

32 

13 

41 

36 

24 

12 

61 

45 

29 

1 

! 45 

35 

13 

i 

1 

i 45 

36 

20 

13  i 

33 

29 

22 

22 

22 

9 

1 29 

27 

19 

14 

38 

26 

21 

26 

19 

9 

26 

26 

21 

15 

20 

18 

14 

1 

1 

1 12 

13 

5 

1 

18 

17 

14 

16  i 

1 16 

14 

11 

1 

; 9 

10 

3 

[ 14 

13 

11 

17  1 

1 30 

28 

21 

! 20 

21 

8 

28 

1 

27 

0 

18 

54 

40 

i 

' 38 

43 

20 

1 

' 54 

55 

0 

.?■ 
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Table  11 


5 

6 

7 

8 

23 

17 

12 

11 

12 

6 

1 

0 

18 

18 

0 

0 

28 

22 

14 

13 

16 

9 

2 

1 

24 

22 

0 

0 

18 

15 

12 

11 

9 

4 

1 

0 

19 

0 

0 

0 

20 

16 

13 

12 

10 

5 

2 

0 

21 

0 

0 

0 

10 

11 

9 

9 

2 

2 

0 

0 

0 

0 

0 

0 

8 

9 

8 

7 

1 

1 

0 

0 

0 

0 

0 

0 

18 

21 

17 

15 

6 

7 

3 

1 

0 

0 

0 

0 

35 

39 

31 

28 

15 

17 

9 

6 

0 

0 

0 

0 
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1.  Program  NAVYS  (FORTRAN) 


I Jd  I '(  -j,  1 >,  5 6,  ■? ) , ^ 7 ) , -^LC  6, 7 ) , 7 ) 

1 do  Li:'EdSIJ'J  I^L(  o,  7 ) , I SL(  1 d , ‘4.' , F3;(  a)  , F SLC /4) 

1 1 U CALL  J'^E.'JFC  1 , "^ILl  " ) 


Ud 

CALL  d'’:;  JF  ( d,  "flLJ  ■■  ) 

1 jU 

CALL  d^C'.'F  C d,  "=>IL4‘' 

UU 

CALL  jP£:dF(-4,  ■•’ILo"  ) 

1 o 0 

CALL  0->£;\IF(o,  "PIdd") 

■ 

1 c d 

CALL  J^C.dFCo,  ■’’I  30") 

17d 

7= 1 . 7o 

1 3d 

o A = . J 7 3 

1 dd 

■;2  = . Udo 

£iUU 

0 = 0 

21  J 

:<=l  i 

-.-id 

<;=  5 

;:.3  J 

DO  1 1=1,  K 

LAd 

P£AD(l,ld>  ( I <(  1 ,0,L),L=  1 , <> 

iioU 

'^CADC  d,  1 d ) ( I <C  d,  J,  L)  ,L=  1 , 

h 

Cod 

OEADC  d,  Id)  C I ■:<  d,  J,  L)  , L=  1 ,M) 

■ 

c;7  J 

PZADC  4,10)  C I ■:(  4,  J,L)  , L=  1 ,--l) 

• 

diJ 

'3£AD(J,1U)  ( I '(  0,  J,  L)  ,L=  1 , M) 

J ^ (j 

=.£AC(o,  Id)  ( I<(u,  J,L)  ,L=  1 ,:3) 

jjd 

1 J 

FO.^.-iATC  ,3  13) 

old 

1 

CO  JTI.JUZ 

32U 

DO  1 7 0 0=1,4 

do  J 

DEADC  3d , 1 7o)  ( I ?L(  I , 0) , I = 1 , 1 

o-^d 

1 7c 

FOliO  ATC  1 U I 1 ) 

did 

173 

CO  JTI  10£ 

• 

1 

do  d 

,10  = 0-  1 

J?  J 

LO  177  1=1,4 

J3  J 

F 3 'K  I ) = .). 

J5»  J 

F.'LC  I ) = d . 

-*0  J 

1 7 7 

CO  J no  J£ 

■ 

44  1 d 

DO  1 U 0 1 = 1,  M 

~tdJ 

I T£C -30,  1 1 II ) I 

■40  d 

1 1 J 

FO?:iAT<  d1,  I 3 ) 

1- 

j 

DO  o 0=7,  < 

■ 

‘4  3 d 

Dd  7 L=l,.11 

d 

L L = L ♦ 1 

^7  d 

CO  TO  dJ-J 

^ jU 

C.U  U 

0=  J . 

49  U 

V=  J . 

o’OO 

DO  3 0C=1,'< 

J 1 d 

IF< I 1C  I , JS.LS)  ) 3,3,4 

n 

320 

4 

:j=-jtl  <(  I , 05,  LS) 

I 

odO 

V=U^I <( I ,OS,LL) 

1 

• 

- 22  - 

4 

.4 


0 -4  J 

J 

coJTi  J J.-; 

5 IJ 

•"  ( I , L ; = c •;  * . 3 > ' c :j  ♦ 1 . ) 

o7  G 

v=33-3r  ( C '’C  I , L)  / ( 1 . - PC  I , L>  ) 

53  0 

y=^.  -^ATAMC  R) 

5?0 

S£=5TRT( 1 ./U) 

60G 

CL=Y-7»se 

fc  1 0 

cu=y+?.'»se; 

O ^ J 

'’L  C I , L)  = C : I J ( CL  •'  3 . ) ) ' ‘ i ) 

o 3 j 

■’'.K  I ,L)  = C C 31  JCC  j/  - . ) ' ‘ 

U 4 0 

IFCIFL)  JoJtJ5j,Jod 

Cl  j 0 

Jg  J 

IF  ( r<(  I ,J,LL)  3,5,4 

oo  U 

3 

L 3 “ L 

c7  J 

1 J 

I F ( r<  < I , J , L r>  4 ) 1 i , 1 1 , I c 

o i G 

1 1 

LS=LC-  1 

l;  4 J 

o3  T4  1j 

7 G G 

i - 

T<-IY(  T , J,L'3> 

/ J D 

'J'J=T  C 

7 1 J 

LF^L-  1 

7 ^ J 

IFL=  d 

7 oG 

G J T J .iJJ 

7s.  J 

3 o G 

JC  I ,L) 

; J J 

rL  = =>L,C  I , L) 

V -J  G 

I^'IC  I ,L)=K3:JCG.',,.3U,  TU.' 

/ / J 

I DLC  I , L>  = < iLC  G , ’ i,  TL> 

/3g 

G J T J 7 

/ 4 J 

4 

L5 - L 

7 4 3 

J-t=  I CC  I ,.J,LS; 

i J'J 

I FL=  1 

> 1 J 

GJ  T3  -'Go 

> .7  il 

JO  U 

TU=-=‘JC  I , L) 

.)  J J 

TL=^LC I , L ) 

3 4 J 

ip-.'c  I ,i-)  = :'C3'jc  (.A,  ;'JJ,  t:;4 

i 30 

I^LC  I , L)  = r<3L(  GR,  J:J,  TL  4 

3 3 G 

7 

CJ.NiTIJL'G 

i 7 d 

VRIT-C  ou,  J J ) I,  ( I •’■JC  I ,L)  , L - 1 , < 3 

3 3G 

JG 

FO  R.IATC  y<,  I 3,  7 I j4 

IVG 

WRITCC  oo,  0 1 4 ( I OL(  I ,L4  ,L=  1 , -l-n 

4 jO 

J1 

FJR'1  MC  <<,  7 1 j4 

4 1 G 

•;RI  TCC  5o,  3 1 4 C I Y(  I , J,  L4  , L=  J,  14 

4-0 

LJ  oGo  IL=l,4 

4JG 

LF=I3L<J-  i,lL) 

4<«G 

1FCLF4  oGG,oGG,ol 3 

4 j d 

c 1 d 

FS'JC  ru4  = FSU(  IL44  Ioij(  i ,lF4 

4 b G 

F3LC  IL4=FFL<  IL4-U'3L(  I ,LF4 

47  G 

G J 

CJ.nJTI  \I'J£ 

4 3G 

n 

CJ'-iTIJ-JF, 

4 4U 

1 Gd 

CJdTi  J';e 

1 U J 0 

VR I T C ( o o , 7 G J 4 

1 G 1 G 

70  J 

FOR  IATC  1 d <,  1 i ITJTAL3  ='J-ECA3TI 

1 UJJ 

'•'PI  TLCoo,  7 1 U 4 C J-  R'JC  I 4 , 1=  1 , 44 

1 UJO 

VR  I To  ( oo  , 7 1 J 4 C F RL  C I ' , I = 1 * 't  4 

1 (J44J 

7 1 J F ) ^-IaTC  0 <,  '^Fl  J . t ) 
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1 U Jii 

•:  JO 

1 OoO 

.-'J.-Ji:TI  3 1 K.iLC  M,  T.' 

1 a?  j 

r,-  n 

A=  1 

1 U-JU 

u-  T 

1 1 Uu 

J=0 

1 1 1 J 

1 

.-■J=C  U.J,r-I,  V) 

i 14J 

IF<F.J-C'.)  3,2,  < 

1 1 jU 

J 

J = J*  1 

1 i u u 

GO  TJ  1 

11  7 0 

ci 

I F C J - 1 J 

1 1 

j 

J = J-  1 

1 1 j J 

o 

K3L  = J 

1 <iUU 

^£T;J3:-J 

1 dl  U 

£ND 

1 

fuaictiq;\i  =,  vi,  t) 

1<£3U 

G=P 

1^40 

12o0 

V*T 

I2o0 

J*  0 

1270 

1 

FJ  = C3( J,M,  W) 

1300 

IF(F.I-G)  2,3,3 

13  10 

ri 

J = .J+  1 

1 320 

GO  TO  1 

1 330 

J 

KB'J=0 

1 3^0 

RGTUR'l 

1 3 jO 

e:jd 

1 360 

rjMCTIO'J  C3(I,\I,^) 

1370 

0=  I 

1 330 

II 

1 3?0 

j=D 

I 400 

IF<.<-  3 j)  11,11,12 

l4l  U 

1 1 

AM  = ^ 

1 420 

3= 1 . - 0 

1430 

I F ( J ) 1,4,3 

144  0 

1 

c:3=o. 

1 4 j 0 

GO  TO  10 

1 4o  0 

<L 

C3=P*».^ 

147  U 

GO  TO  10 

14^0 

J 

3=3‘»»  '1 

1430 

C 3=  3 

1500 

00  4 -(=  1 , J 

1 3 1 0 

AK=< 

1 520 

A1-I-4K+  1 . ) / C R*4iO 

1 330 

C3  = CB^  3 

1 340 

COWTI JUE 

1 550 

GO  TO  10 

1 3b  U 

1 d 

II 

1 570 

(VJ=.; 

1 j30 

’:=  C A.J+.  3- AM  "O)/  SORTC  AM>  ( 1 . 

24 
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I :>-JJ  i:  JiK(  ^ ' 

loJll  I ii  I'cf  ''(J 
lulu  ^ J u 

lic’O  MjiJJTl  J..'  CJL-:(V> 

l3oU  <=Y 

I3^U  ISVTCI=0 

l3oJ  IrCO  I , £l,  d 

I3t)j  I <=<\BSC<) 

137  0 I?VTC3= I ' 

1330  ^ ==  . I 1 

13J3  31=.3lij3lD3 

1300  3<i=-.35c5o373 

13  10  oo=  1 . 73  1 -47  7 3 

13^0  .3<4  = -1.3eld553 

1330  35=l»j3  0«i7‘4*i 

13^0  T=  1 . /<  1 . +'3‘*'0 

1 330  . 3 3 33<4<i3*£X'3( -X^  </ =i.  ) 

1 36  0 3 3D<=  1 . - ’3*  ( -,1  32*T*T>  33'*  < *3  ) ♦H4'*  ( T*  »^)  + 35*  ( T**5  ) ) 

1 370  IFCISVTCi;  3,3,3 

1330  3 0MD<=  1 . -QfJDX 

1 330  3 C:'JDX='3MI)X 

cOUO  Ti'TLJ'lN 

<d  0 1 0 £ J D 
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2 • Program  NAVYR  (FORTRAN) 


1 u .)  M I -I  - I J i ! < ( - , 1 5 , ^ C , / , 7 > , DL  C 6 , 7 ) , I D'j  < o , 7 ) 

1 J o i.I  '!!•:  vJS  I J.O  I ■’uC  , 7 ; , I ?L(  1 ij  , ‘-i  > , ? =:'J<  4 ) , r 3L(  ^ ) 

i I u c M.L  3'’;-: 'IF  ( 1 , 3IL 1 ■■ ) 


1 'J 

CALL  r=Z  JF  ( ^ , ”ni  L J ■ 

A 

1 J J 

Lv'aLL  C ji  ’^IL'L 

> 

CALL  3^1  J.- ( -i,  lL'-  ' 

> 

1 J u 

C-'.LL  j j ■ 

> 

1 u J 

CALL  J^C.-'IFC  o,  "'31  J v’ 

) 

1 7 J 

^=\  . 7 0 

1 <0 

J - '0 

1 7 J 

1 J 

^ JJ 

•:=  J 

^1  0 

L J 1 ■ J - i , K 

' 'j 

'’t'.ACC  1 , 1 J J ( I •(  i , 1 , 

L) , 1 , 

•1 ) 

i J J 

•3£mC(  , 1 0 > ( I 

U ' 1 , 

1 ' 

il 

'■■'  c.  A J C .3 » 1 ij  > C I - 1 . 1 f 

L)  1 , 

1 ) 

^ J J 

J , L = 1 , 

<i'jO 

Ri'AuC  0,  1 J)  ( I <(  0,  J, 

L)  ,L=  1 , 

1 ) 

/ u 

■3  CALC  n , 1 0 > ( I < C , .f , 

L) , L=  1 , 

1 ) 

- i J 

J o 

F JR  < AT  1 J I ->  ) 

1 

C J.'JT  I J'J£ 

iJJ  17j  -1=1, 

/ 

"f  Al3(  j'j,  1 7o)  C I AL(  I 

,17,1=1 

. 1 

y i 

i 7 

; FO  iMATC  1 d I 1 

1 / o 

C'J  JTI  J'I£ 

J b J 

i-n-  -1-  1 

J0<i 

D)  177  1=1,. 

J U 

? 5 ■ K I ) = '1 . 

JO. 

F5L(  I ) = ;i. 

J J J 

1 7 7 

CJM'I  J'.'I 

J 1 li 

J D 1 lI  d 1 = 1 , d 

J.' u 

VP.I  Tr.(  00,  1 1 J I 

3 J J 

1 1 Li 

F,)='d.AT(  J<,  I > 

Joo 

DJ  V •)=<,:•< 

J *»  0 

iJO  L=  1,1-1 

Jo  J 

LL  = L*  1 

J 0 0 

C 3 7 J J u 0 

Jo  U 

<:  D J 

J " C*  • 

J70 

V=  J . 

JJU 

oO  J . 1 '••  = 1 , < 

J y u 

IFC  I<(  1 ,.J3,LAJ  J,  J 

, 

4Uu 

H 

U=-j*  I <<  I ,.IS,  L3) 

aJ  J 

77  = V»  I C(  I , Jr.,LLJ 

U 

i 

lC3:-JTI.1'Jc. 

26 
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^ J 0 

^ ( I , 1 . - > ' c ► 1 . 

■*  1 

>'^TC  ( 1 , ' s 1 . - 

MO'J 

v=  , . 1 at4  J ( ■■= .' 

J 

S^=  >7TT(  l ./-I) 

'•*  7 ij 

CL  = Y-7  ► 

^ i J 

C'j=v+7i>s-: 

-*  V 0 

DLC I ,L)  = ( C 51  0(CL/£. 

SOU 

puc  I ,L)  = ( ( srjcc'j/c;. 

o 1 0 

I F ( I FL)  JsUijSO.SoO 

o?  0 

J JU 

IFC  I <(  I , J,LL)  ) 5,5, 

oU  0 

L5  = L 

o 1 J 

i j 

I r C I < C I , J , L 5 i 1 , 1 

C.  Cl 

1 1 

L " = L ~ - 1 

JJU 

.',J  T'J  1 J 

U iJ 

i .i 

T •=!  <(  I ,.I,LS) 

6 ^ 0 

L5=L-  1 

w/  J J 

I FL=  J 

u j7 

iJ  TJ  ^-')J 

OO  II 

: J J 

j-3=T<^'>-!C  I ,L) 

c)  7 J 

T'jL-  1 . -'’U<  I , L) 

J 5 li 

7P='ji?  "piJL 

,>U'J 

"'’=3  7 5T(  V^’) 

7 JO 

■•r>+  * 3"’ 

7 1 U 

I ( T'<-V  = 5 

7d0 

1 

I P J ( I , L)  = I .J  T C T 

7o0 

(;j  T) 

i’  0 

I ■•■’(  I 

7j0 

- J 

•JP=T  'PLC  I , L.' 

7u0 

?LL=  1 . - •’LC  I ,L) 

7 7 J 

Vr>=  'jn  ‘n^L 

/i  0 

S'’=$r'.T(  v">.' 

/ J 

v=  IJ")-  3 ■’ 

ill  J 

IF(T<-UP;  c-i,.ij,Yj 

il  0 

I ’’L  C I , L ) = I !J  T ( T ) 

SdO 

■)J  T.j  7 

i J u 

>0 

i^Lc  I ,L)  = 

3s.  J 

5J  TJ  7 

i J J 

) 

L3  = L 

5<:U 

T C=  l YC I , J, LS7 

5 7 J 

IFU=  1 

5 J J 

5 5 TJ  iiijb 

J J 

3 J li 

"'■JL  = PIJC  I , L 

U 1 0 

•JJ=-r<  V JIJL 

Uc:  U 

VP=|JD  ► ( 1 . - V!L  > 

UiJ 

3n;=  7Q> 

V<4  il 

i/p  = :,"■)♦  );  + 3 p 

7j  J 

I F ( T<-  .;p;  J 1 , 0 1 , J . 

>o0 

-5  i 

I'  j<  I , I.  > = 1 JT(  T C) 

J7U 

1. ; J T J J 5 

7 50 

J i 

IPJ<  I,L7  = i:OT<'’P7 

27 


T-358 


^ i ,) 

3 i 

■L'.-  •’Lv  I , w 

1 J 0 J 

• J ~)  3 T ' ^ " L.  L 

1 0 4 J 

V'=  U 1 . - (LL) 

J J j 0 

5 3^3  ^ " r ( V ^ ■' 

1 UoD 

;jr>_  7 >^30 

1 0 / u 

I F cr / 1 , 7 i 

f 

1 j Ui 

! 1 

I ^Lc  t , i. ' ^ I jrc  r < 

1 u y j 

w J T j 7 

1 1 b u 

/ 

I^l.C  I ,L>  = I JTC'--^ 

A 

1 1 1 J 

/ 

c:j  JT I riz. 

1 1 y ij 

••■^I -Tic  ,.o  , J 1 7 I, 

( DO’C  I , D , - = 1 , 1 ^ ? 

1 1 OC 

J 

Ii  F J {.-l  A'C  J<,  1 3,  / 

1 J ' 

1 1 .tU 

WR  I TcC  , J 1 ) (I 

•^LC  I , D ,L=  1 , ‘1  •) 

1 1 Ij 

> 

1 .-  y- 'i.xrc  5 <,  7 I ) ) 

1 1 >IJ 

I J 1 ) Cl 

<(  I , I , , L = . , -'  ? 

1 lu.-: 

. ) J 0 J 1 L= 1 , » 

1 1 ..  J 

LK  = I SLCJ  - ' , ID 

M n4 

1 !■  ( L.Z  ) (1  J il  1 0 7 0 1 

3 1 ij 

i 1 J 

1 iJ 

7MC  D7  = -F';C  ID 

+ D'K  I , OK  ? 

1 1 - 3 

F ?,LC  ILDFSLC  IL? 

* I ">L  C I , LF  ? 

i U:  i 

.J 

:j  jti  j7j : 

1 1 7 J 

3 

C..;  JTi  rjL 

) J i J 

1 4 J 

LJJTI  JUZ 

i 1 1 J 

'..'^ITLCoc  , 7b0) 

J 1 io 

7 0 J 

FUl-lATC  1 JX,  1 iilTDTALS  F J ’•■ICAS’' I IG? 

1 1 <7 

' ’ I T £ ( 0 0 , 7 1 J ) 

( FSUC I ) , 1= 1 , 4? 

1 1 i i 

-7  I T £ C 0 0 , 7 1 u ? 

C FFLC I ) , 1 , 4? 

1 1 i 7 

7 1 

J F 1 ATC  3 : , -4h  1 J 

. b ? 

1 ) y J 

•:  JO 

28 


r-'X  ' 


THE  GEORGE  WASHINGTON  UNIVERSI  IT 
Program  in  Logistics 
Distribution  List  for  Technical  Papers 


The  George  \\ashing(on  Universily 
Office  of  Sponsored  Resetrch 
Library 

Vice  President  H.  K.  Bright 
Dean  Harold  Liehowitz 
Mr  J hrank  DouMeday 

ONH 

Chief  of  Naval  Research 

(('odes  200.  430D.  I021P) 
Resident  Representative 

OPNAV 

OP  40 

D('N(>,  1 ogistics 
Navy  Dept  Library 
OP  *^1  I 
OP  <*64 


Army  Logistics  Mgmt  ('enter 
Port  Lee 

Commanding  Officer.  USALDSRA 

New  Cumberland  Army  Depot 

US  Army  Inventory  Res  Ofc 
Philadelphia 

HO.  US  Air  Force 
A FADS  3 

Griffiss  Air  horce  Base 

Reliability  Analysis  Center 

Maxwell  Air  Force  Base  Library 

Wright-Patterson  Air  Force  Base 
HQ.  AF  Log  Command 
Research  Sch  Log 


Naval  Aviation  Integrated  Log  Support 
NAVCOSSACT 

Naval  Cmd  Sys  Sup  Activity  Tech  Library 

Naval  Flectronics  Lab  Library 

Naval  Facilities  FngCmd  Tech  Library 

Naval  Ordnance  Station 
Louisville,  Ky. 

Indian  Head.  Md. 

Naval  Ordnance  Sys  Cmd  Library 

Naval  Research  Branch  Office 
Boston 
('hicago 
New  Yt>rk 
Pasadena 
San  Prancisco 

Naval  Research  I ab 
Tech  Info  Div 

Library.  Code  2029  (ONRL) 

Naval  Ship  Plngng  Center 
Philadelphia.  Pa. 

Hyattsville.  Md. 

Naval  Ship  Res  & Dev  ('enter 


Defense  Documentation  Center 

National  Academy  of  Science 

Maritime  Transportation  Res  Board  Library 

National  Bureau  of  Standards 
Dr  E.  W.  Cannon 
Dr  Joan  Rosenblatt 

National  Science  Foundation 

National  Security  Agency 

WSKG 

British  Navy  Staff 

Logistics,  OR  Analysis  Establishment 

National  Defense  Hdqtrs,  Ottawa 

American  Power  Jet  Co 

George  Chernowit/ 

ARCON  Corp 

General  Dynamics.  Pomona 

General  Research  Corp 
Dr  Hugh  Cole 
Library 

Planning  Research  ('orp 
Los  Angeles 


Naval  Sea  Systems  Command 
Tech  Library 
Code  073 

Naval  Supply  Systems  ('ommand 
Library 

Cap!  W.  1 . Nash 

Naval  War  College  Library 
Newpj»rt 

BUPFRS  Tech  Library 
FMSO 

Integrated  Sea  Lift  Study 
USN  Ammo  Depot  Earle 

USN  Postgrad  School  Monterey 
Library 

Dr.  Jack  R.  Borsting 
Pr<»f  ('.  R.  Jones 

US  Marine  Corps 

('ommandant 

Deputy  ('hief  of  Staff,  R&D 

Marine  Corps  School  (^antico 
Landing  Lorce  Dev  Ctr 
Logistics  OfHcer 

Armed  l-orces  Industrial  College 
Armed  Porces  Staff  College 

Army  War  ('ollege  Library 
('arlisle  Barracks 

Army  (’md  & (ien  Staff  College 

US  Army  HQ 

LTC  (»eorge  L.  Slyman 
Army  Irens  Mat  Command 


Rand  Corporation 
Library 

Carnegie  Mellon  Cniversily 
Dean  H.  A.  Simon 
Prof  (>.  I'hompson 

Case  Western  Reserve  University 
Prof  B.  V.  Dean 
Prof  John  R.  Isbell 
Prof  M.  Mesarovic 
Prof  S.  ^acks 

Cornell  University 

Prof  R.  E.  Bechhofer 
Prof  R.  W.  Conway 
Prof  J.  Kiefer 
Prof  Andrew  Schultz,  Jr. 

Cowles  Foundation  for  Research 
Library 

Prof  Herbert  Scarf 
Prof  Martin  Shubik 

Florida  State  University 
Prof  R.  A.  Bradley 

Harvard  University 

Prof  K.  J.  Arrow 
Prof  W,  (J.  Cochran 
Prof  Arthur  Schleifer,  Jr. 

New  York  University 

Prof  ().  Morgenstern 

Princeton  University 

Prof  A.  W.  Tucker 
Prof  J.  W.  Tukey 
Prof  Geoffrey  S.  Watson 


NATIONAL  ADVrS< 
NATIONAL  BOftEAl 
AMERICAN  SOCjr 


or  crvii 


ameri^n 

THE  AMERlCAN^eRgreTTjSSti 
'-THE  SOCIETy  OFJN^MSyir. 
AMERICAN 

OI STRICT  or  COLlffBTX  «)Cf 

CIFE  C3rC»IMU£ 

wKSftm  ^ 


rayiEE 

seiEN 


:l2.U:i^“\^rU.'K.CLL 


UlUiLVV.  aM>UL>i 


Wl-UW-UlsUY 


'/W'12l«UVW3VU;‘ 


.rvnwi 


I ( ) ( ( i|)(’  u III)  llic  cxp.indiii!’,  I('(  litK )l( n;\ , nm  si  k m>I\  imisl 
he  .issuM'd  III  ,1  ( oiiliMuin.i;  suppK  nl  ni;( )ii ms|\  li.iini'd 
,ind  I'diK.ilcd  ('M^mci’ts.  Ilir  S(li()iil  ol  I mpiu'iMin.v;  ,ind 
Applx'd  SiiciKc  IS  ( ( jinpli 'll ‘l\  I immiincd  In  ihis  nh 
|('(  llU', 


